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. Abstract 

We evaluate the discovery reach of the ATLAS experiment for down type isosinglet quarks, D, using both their 
neutral and charged decay channels, namely the process pp — * DD + X with subsequent decays resulting in 2£ + 
2j + rfj, 3£ + 2j + If j and 2i + 4/ final states. The integrated luminosity required for observation of a heavy quark 
is estimated for a mass range between 600 and 1000 GeV using the combination of results from different search 
channels. 



1 Introduction 



In an effort to unify the forces known to the Standard Model (SM) within a Grand Unified Theory (GUT), one is led to 
consider extended families of fermions. One of the often quoted GUT symmetry groups having SUc(3) x St%(2) x 
Qh' Uy(l) of the SM as a subgroup is the exceptional Lie group with its 27 dimensional fundamental representation 

Q],ll|. Such a large symmetry group not only puts the known quarks and leptons in the same family, but also predicts 
additional particles, such as iso-singlet quarks and new heavy gauge bosons. In the literature, the new charge q=— 1/3 
iso-singlet quarks are denoted by the letters D, S, and B, representing the additions to the first, second and third SM 
families, respectively. 

Here, we investigate the discovery potential of the ATLAS experiment for the down type isosinglet quark of the 
first SM family (D quark) at the forthcoming LHC accelerator. In accordance with the observed quark mass hierarchy, 
we will assume the first family's isosinglet quark D to be the lightest of the particles predicted by the £g group and 
hence the first one accessible at the LHC energy. It is assumed that the other new particles predicted by the model, 
being very massive, do not contribute to the D quark search. The present experimental limit on the mass of such an 
iso-singlet quark is mo > 199 GeV [3 1. A second assumption will be that the intra-family mixing is stronger than the 
mixing between families. Although we will concentrate only on the D quark, the results will stay valid also for the S 
quark, as jets from first and second quark generations are experimentally indistinguishable. The investigation for the 
B quark requires detailed information about the detector response to b-jets which is beyond the scope of this work. 
The Lagrangian relevant for the weak interactions d and D quarks is given as J4] : 



sj\%a em 
2V2sin 9w 
\l\na em Tsin^cosi/) 
4sin0iv cos(?w 



u e Ya (1 - Y5)dcos(j) + u e Ya (1 - y 5 )£>sin0 
dy a {\-Y5)D 



W a (1) 



[Dya {4- sin 2 0w- 3 sin 2 <p (1 - y 5 )) D + dy a (4sin 2 0w- 3 cos 2 <p (1 - y 5 )) d] Z a +h.c. 



12cos0M/sin0n/ 



where the superscript 9 represents the usual CKM mixings taken to be in the up sector for simplicity of calculation, 
6w is the weak angle and <p is the mixing angle between the d and D quarks. This mixing is responsible for the decay 
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Figure 1 : Main diagrams for the pair production of the D quark at the LHC. 



of the D quark. The limits on <p can be obtained from the current precision measurements for the 3x3 CKM matrix 
elements assuming that its 3x4 extension has the sum of the squares of the elements of a row equal to one. The 
evaluation of the present values yield | sin 1 < 0.045. Although the upper limit is used throughout this work, the 
results are essentially insensitive to sin as the pair production considered here proceeds mostly via gluon exchange. 

The SM fermions acquire their masses through their interactions with the iso-doublet Higgs field. The Higgs 
mechanism can also be preserved in the £g group structure as an effective theory, although other alternatives such 
as dynamical symmetry breaking are also proposed 0. On the other hand, the origin of the mass of the new quarks 
(D, S, B) should be due to another mechanism since they are isosinglets. For this reason and following the literatureQ, 
Eq. ([T} does not contain Higgs related terms. However, if the mass of the Higgs boson is smaller than the isosinglet 
quark mass and if the interaction between the Higgs field and d quark is considered before the spontaneous symmetry 
breaking but after d — D quark mixing, a new D quark decay channel becomes available: D — > lid. The effect of such 
interaction, depending on the Higgs boson and D quark masses, is to reduce the BR(Z) — > W u) from 66% down to 
50% and BR(D -> Zd) from 33% down to 25%, giving way to BR(D -> hd)=25% (for details see @ and Q). 



2 Analysis 

The leading Feynman diagrams for the pair production of the D quark at the LHC are given in Fig.Q] The Lagrangian 
in Eq. (HJ was implemented in two tree level generators: CompHep (8) and MadGraph 0. The pair production signal 
cross section was calculated in both generators as a function of D quark mass, yielding practically equal cross sections 
for the same generator level cuts. The result from CompHep is plotted as a function of the new quark mass in Fig. |2] 
Contributions from the electroweak interactions to the production cross section are found to be much smaller than those 
from the strong interactions. Below a D quark mass of about 700 GeV, the dominant contribution to the production 
cross section is from the diagrams in the top row of Fig. |2]but as the mass increases, the quark fusion contribution 
takes over according to CTEQ6L1 PDFs [ 1 1 . One should note that single production of D quarks (in association with 
a light SM quark) is also possible in this model. However, since this type of subprocesses, in comparison to the pair 
production, would give a different final state the production of single D quarks is considered in a separate study ifTTI . 

The £g GUT model does not predict the masses of the isosinglet quarks. Therefore, this study scans some plausible 
values for the D quark mass (600 - 1000 GeV) to estimate the experimental reach of the ATLAS detector lfT2l . With 
such a large mass, the isosinglet quarks are expected to immediately decay into SM particles. A number of decay 
channels of the D quark pairs are summarized in Table Q] The process in the first row has been previously studied 
and found to be feasible for discovery up to a mass of about 1 TeV, for an integrated luminosity of 200 fb -1 lfL3l . 
The present work aims to evaluate the possibility of discovery of the D quark in the production channels shown in 
the remaining rows of the same table using the final states in column 3. The high transverse momentum of the jets 
coming from the D quark decays can be used to identify the signal events. As for the backgrounds, all the SM tree 
level processes yielding two jets and two vector bosons (Z Z and W — Z) are considered. In this feasibility study, the 
detector related background events originating from jet combinatorics and misidentification such as two vector bosons 
plus one or three jets are not taken into account. The background cross section was calculated only in MadGraph as 
this generator is suitable to be used on large computing farms. In the numerical calculations, we set the QCD scale of 
the signal to the mass of the new quark while keeping it as the mass of the Z boson for the background, yielding more 
conservative results compared to the investigation in lfl3l . For example, the 3<7 observation limit gets reduced by 15% 
in the in the Al + 2j channel studied in that work. The simple requirements imposed at the generator level are: 

|»lp I < 2.5, (2) 
PT, P > 10 GeV, 
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Figure 2: Tree level pair production of D quarks at LHC as a function of the new quark's mass. 



Table 1: List of the studied signal channels: the last column contains the total branching ratios to the final state 
particles where I represents ji or e . 



DD^ 


Final State 


Expected Signal 


Decay B.R. 


Total B.R. 


ZZdd 


z^iiz^il 


4 £+2 jet 


0.07x0.07 


0.0005 


0.33x0.33 


z^eiz^vv 


2£ +2jet + E/ T 


2x0.07x 0.2 


0.0028 




Z->etZ^ qq 


2£ +4 jet 


2x0.07x0.7 


0.0107 


ZWdu 


Z^££W ->/v 


3£+2 jet+E/ T 


0.07x0.21 


0.0065 


2x0.66x0.33 


Z^UW -^qq 


2£+4 jet 


0.07x0.68 


0.0211 



R P p > 0.4, 
\r]z.w\ < 5.0 

where R is the cone separation angle between two leading partons (p) giving rise to two jets. Selection of the pseu- 
dorapidity region is driven by partonic spectra distributions, which are peaked in the central region. Decays of Z 
and W bosons were performed by PYTHIA fBl which handled the initial and final state radiation together with the 
hadronization. The detector response was estimated using the ATLAS fast simulation program ATLFAST lfT31l within 
the offline analysis framework, ATHENA vl 1.0.41 . 

2.1 Search using the 2£ + 2j + $t channel 

This search aims to exploit the relatively large branching ratio seen in second line of Table [T] However the invisible 
decay of one Z boson, makes reconstruction of both D quarks impossible. The expected viable signature is two 
isolated leptons, two energetic jets and large missing transverse energy. The reconstructed invariant mass from the two 
leptons is required to yield the Z boson mass with an error of ±20 Ge V which is an order of magnitude larger than the 
reconstruction resolution ±2.5 GeV. The two most energetic jets are selected as candidates for the parton daughters 
of the D quark decays. However, the ambiguity in the association of a jet and a Z boson, makes each event contribute 
twice to the D quark invariant mass histogram, one with the correct association and one with the wrong one. The 
background for these events is estimated by taking into account all the SM processes which give two energetic jets 
and two Z bosons, allowing one of them to decay invisibly and the 'other one into electrons or muons. The total cross 
section of such processes with the generator level cuts in Eq. except pj p > 100 GeV is 2.65 pb. The other SM 
processes with 2£ + 2j + background are expected to be negligible, given the high Ifj in the topology of the signal 
events, and the Z boson mass reconstruction requirement. 

The simple cuts for this channel are presented in the set below: 

#(e,n) = 2 (3) 
E/ T > 150GeV 
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Figure 3: Invariant mass spectra of signal and backgrounds for different values of D quark mass in 21 + 2j + Fjj 
channel 



Table 2: The expected number of events in the mass window around the peak and the statistical significance for 
2£ + 2j + E/t channel after the event selection cuts in equation set[3jor an integrated luminosity of 100 fb~ 1 . 



Mass (GeV) 


600 


800 


1000 


S 


53 


19 


4 


B 


12 


13 


5 


significance 


15.3 


5.3 


2.1 



M a = 90 ±20 GeV 
#jets > 2 
Pt jet > 150 GeV 

The resulting invariant mass plots for different D quark mass values are given in Fig. [3] for one year of LHC 
data taking which corresponds to an integrated luminosity of 100 fb . The solid histograms represent the events 
originating from the SM background and the dashed ones from the signal. Both distributions were fitted separately to 
reduce statistical fluctuations and to better estimate the total number of events in the histograms. The number of events 
for both signal and background were computed by integrating the fit functions in the range of ±2<7 from the center 
of the Gaussian fit to the signal peak. Table [2] contains the number of signal (S) and background (B) events together 
with the signal significance calculated in this way for an integrated luminosity of 100 fb -1 . In this note, the signal 
significance is calculated as S/yB if S + B > 25 , otherwise Poisson statistics are used to obtain the probability P, of 
compatibility with the background. In the latter case the significance of the Gaussian distribution yielding the same 
probability, P, is given. 

2.2 Search using the 3£ + 2j ' + $j channel 

This channel has one D quark decaying through a Z boson and another D quark decaying via a W boson, i.e. DD — > 
ZW jj final state. The Z boson is expected to be entirely reconstructed from two electrons or two muons. The W boson 
can also be partially reconstructed if a good measurement of missing Ej can be achieved. The signal is sought in the 
final state made of a pair of energetic leptons (electrons or muons), accompanied by a third high pj lepton (muon or 
electron), two high pj jets and the missing transverse energy. 

For the background studies, events consisting of the SM WZ + 2j processes were generated. The background 
cross section before event selection is 2.30 pb from the W~ case and 3.89 pb from W + case, summing up to 6.19 pb 
with the generator level cuts in Eq. $2& except pj. p > 50 GeV. The composition of background final state and allowed 
decay modes of Z and W bosons were chosen to be exactly the same as for the signal events. The other possible 
background contributions due to misidentified or undetected leptons, such as ZZ + 2j or WW + 2j events, will be 
significantly suppressed due to the event selection: rather high transverse momentum cuts on jets and the requirement 
of reconstructing both Z and W bosons in the same event. 

Two cases were separately considered during the event selection for the lepton composition of the final state. These 
two cases and the event selection cuts for a D quark mass of 600 GeV are: 

A. the Z and W decays involving leptons of the same generation (electrons or muons) were selected using 

#e = 3 or #n=3 (4) 

B. the Z and W decays involving leptons of different generation (electrons or muons) were selected using 
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Figure 4: Jet (left), electron (center) and muon (right) transverse momentum spectra for 3£ + 2j + Jf.j channel. The 
shaded area represents the SM background contribution. 

#e = 2 = 1 or #e = 1 #/x = 2 . (5) 

The second case aims to ease the identification of leptons originating from Z and W bosons. The other analysis 
cuts are common to both cases: 



p T j > 80GeV (6) 
p T , e ,n > 20GeV, 
|T?C/,eorju)| < 2.5, 

p'f ss > 30GeV, 
M\f = 90±20GeV, 
M^ c = 80±20GeV. 

For higher D quark mass values, the jet transverse momentum cuts were slightly increased to obtain a favorable 
significance value wherever it was statistically possible in both cases. 

Transverse momentum spectra of jets, electrons and muons in the signal events for various D quark mass in 
comparison with the corresponding spectra from the background events, are presented in Fig. @] after the selection cuts 
defined above. One can see that the signal events feature more energetic jets and leptons than the background data. 
Pseudorapidity distributions of jets and electrons in the same events are shown on Fig. [5] Partons from signal events 
are produced mostly in the central region, more so as the D quark mass increases. Missing transverse momentum 
distribution, produced in the event as a result of the semi-leptonic W decay is shown on the Fig. [6] As it can be seen, 
signal events produce more energetic W bosons, which results in larger missing transverse momentum vector with 
respect to the SM WZjj events. 

The reconstructed invariant mass of the two leptons was required to be in the (90 ±20) GeV mass window to ensure 
their origin from the Z boson. The W boson "visible mass" was reconstructed from the measured lepton momenta and 
the Tf.j in the transverse plane, assuming negligible neutrino longitudinal momentum. A mass window constraint of 
(80 ± 20) GeV was applied to the reconstructed M£" : to select events with a W boson. 

In case A, the identification of the two leptons originating from the Z boson decay forced the use of the third 
lepton to reconstruct the W boson invariant mass. The method used to identify these two leptons and thus to avoid the 
problems due combinatorics related to the issue of lepton assignment to Z or W, was to select those giving the most 
accurate Z boson mass. In case B, as the leptons originating from Z boson were of a different SM family as compared 
to the one from the W decay, reconstruction of the W invariant mass was straightforward. 

As QCD radiation from energetic partons may result in an increased number of produced jets, a correct assignment 
of selected jets might not always be possible due to combinatorics. Moreover, the signal events produce a significant 
number of jets, as can be seen from Fig. [7] where the jet (pr> 70 GeV) multiplicity distributions for the signal events 
of a D quark mass 800 GeV and for the background events are presented. Such a large jet multiplicity could pose 
a problem for the identification of the jet originating from the D quark decay and hence for the reconstruction of its 
invariant mass. In both cases the invariant mass of the D quark, decaying to a Z boson and a jet was reconstructed 
by combining those two leptons with the most energetic (or the next to most energetic) jets only. The ambiguity in 
assigning the two most energetic jets to the two reconstructed bosons can be solved by calculating the invariant mass of 
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Figure 5: Jet (left) and electron (right) pseudorapidity distributions for different D quark masses using the 3£+2j +Efr 
channel. The shaded area represents the SM background contribution. 
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Figure 6: Missing pj distributions for different D quark masses using the 3£ + 2j + $r channel. The shaded area 
represents the SM background contribution. 




Figure 7: Number of jets with pj > 70 GeV for signal (white) and background (shaded) events in the 3£ + 2j + $t 
channel for an example D quark mass of 800 GeV. 
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Figure 8: Invariant mass reconstruction spectra of signal (shaded) and background (hashed) for two values of D quark 
in 31 + 2j ' + tfr channel. The same lepton flavor was found in Z and W decay chains (Case A). 
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Figure 9: Invariant mass reconstruction spectra of signal (shaded) and background (hashed) for three values of D quark 
in 3£ + 2 j +rfr channel. The lepton flavor found in Z and W decay chains was different (Case B). 



the D quarks using all possible combinations and by taking the one which gives the smallest mass difference between 
the D quarks obtained from Z and W in the same event. The D quark invariant mass reconstruction spectra obtained 
with this method are shown in Figj8]for case A, where Z and W decay modes involve leptons of the same generation 
and in Fig|9]for case B, where Z and W decay modes involve leptons of different generations. The figures contain, for 
an integrated luminosity of 100 fb~ 1 the reconstructed invariant mass distributions oiD^Zj and D — > Wj processes. 
The latter have wider invariant mass spectra due to the use of the "visible mass" of the W in the reconstruction routine. 

The number of signal (S) and background (B) events together with calculated statistical significance, corresponding 
to 100 fb _1 integrated luminosity are presented in tables [3] and |4] for type A and B events. The number of events in 
the signal and background processes were calculated in sliding mass windows of ±2 x 8M around the peaks of the 
signal events for both Wj and Zj invariant mass histograms. This selection aimed to select the events where both 
D quarks were successfully reconstructed. The example of the sliding windows mass cuts to calculate the number of 
events in signal and the background is presented in Fig. [I0]for type B events and for a D quark mass of 600 GeV. 
The selection windows for all D quark mass values are given in tables [3] and |4] for type A and B events. The same 
tables also contain the statistical significance of the signal calculated as S/\/~B if the total number of events were larger 
than 25. Otherwise, the signal observation probability was calculated using Poisson statistics and the significance of a 
normalized Gaussian distribution yielding the same probability was reported. 
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Figure 10: Example of the rectangular selection accepting events with two correctly reconstructed D quarks. The Zj 
and W j acceptance widths for different mass values are given in tables [3] and |4] Here the selection for type B events 
of mo=600 GeV is shown. 



Table 3: The expected number of events in mass windows of width 8M around the peak and the statistical significance 
for the 3£ + 2j+$T channel, after the cuts given in[6] for type A (same lepton generation) using 100 fb _1 of integrated 
luminosity. 



Mass (GeV) 


600 


800 


1000 


8M (D^>Z jet) 


120 GeV 


140 GeV 


160 GeV 


<5M (D-> W jet) 


200 GeV 


220 GeV 


240 GeV 


S 


45 


8.3 


0.5 


B 


11.4 


4 


0.25 


significance 


13.3 


3.3 





Table 4: The expected number of events in the mass window around the peak and the statistical significance for the 
3l + 2j + $T channel, after the cuts given in|6] for type B (different lepton generations) using 100fb~ 1 of data. 



Mass (GeV) 


600 


800 


1000 


8M (D^Z jet) 


120 GeV 


140 GeV 


160 GeV 


8M (D— > W jet) 


200 GeV 


220 GeV 


240 GeV 


S 


52 


10 


1.3 


B 


13.5 


5 


0.8 


significance 


14.2 


3.8 
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Figure 1 1 : Multiplicity and pj jets for signal(white) and background(shaded) events for 2£ + 4j channel 



2.3 Search using the 2£ + 4j channel 

This channel has one D quark decaying through a Z boson which is expected to be reconstructed from two leptons 
(electrons and muons) and a second D quark decaying either via a W boson or a second Z boson which can also be 
reconstructed using two energetic jets. One should note that the reconstruction efficiency of a highly boosted vector 
boson decaying hadronically can be lower compared to smaller boosts since the small opening angle between the jets 
can lead to a single jet. The signal events were simulated for both the ZW j j and ZZ jj cases. For this type of process, 
the signal is sought in a final state made of a pair of energetic electrons (or muons) in a multijet environment which 
includes two high p-p jets. The background events consisted of the SM WZ + 2j , ZZ + 2j and Z + 2j processes, 
with generator level cross sections of 6.19, 5.5 and 25 pb respectively with the generator level cuts in Eq. (f2]i except 
p T ,p> 50 GeV. 

The transverse momentum distribution of the leptons in the final state considered in this section do not differ 
much from those presented in the previous section. Therefore, Fig. QT]contains on ly the multiplicity of jets and their 
transverse momentum spectra for signal events using various D quark mass in comparison with the corresponding 
spectra from background events. As it can be seen in the same figure, the signal process produces more energetic jets 
in the final state as compared to the SM background. However the jet multiplicity is higher than the case considered 
in the previous section. In general, high jet multiplicity creates problems in the invariant mass reconstruction due to 
combinatorics. 

The cuts for the reconstruction of the pair produced D quarks with mass of 600 GeV are presented in the set below: 



M{e, n) 




2 


#jet 


> 


4 


PT.j 


> 


80 GeV 


Pti 


> 


20 GeV 


|t?CM)I 


< 


2.5 


M rec 




90 ±20 GeV 


M rec 

jj 




85 ±25 GeV 



For higher D quark mass values of 800 and 1000 GeV, the jet transverse momentum cuts were increased by 10 and 
20 GeV respectively for each next mass value to obtain a better value for signal significance {Sl\/B or from Poisson 
statistics, depending on the available number of events as mentioned before), wherever it was statistically possible. 
The reconstruction of two lepton invariant mass was performed to obtain the mass of Z boson. 

The D quark invariant mass reconstruction spectra with transverse momentum cuts applied are shown in Fig. [T2] 
The mass of D quark, decaying to a Z boson and a jet was reconstructed by combining two highest pj electrons (or 
muons) with a leading, or a next to the leading pj jet choosing the case with the minimum mass difference between 
the two D quark candidates as in the previous subsection. The invariant mass spectra for D — > Z + jet decays are 
shown on the left side of Fig. Q~2] for each considered mass. The decays of a D quark to W + jet were obtained by 
requiring the reconstruction of W mass with two jets in the (80 ± 20) GeV window and then by combination of the 
two "W-tagged" jets with the leading or next to leading pj jet. The invariant mass spectra of D — ► W + jet decays 
are presented on the right side of the Fig. [12lfor each considered mass. The distributions are shown for an integrated 
luminosity of 100 fb _1 . 
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(b) WZjj with m/j=800 GeV 
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(c) ZZjj with m D =600 GeV (d) ZZ jj with m D =800 GeV 

Figure 12: Invariant mass reconstruction spectra of signal (shaded) and backgrounds for 2£ + 4j channel 



Table 5: The expected number of events in the mass window around the peak and the statistical significance for the 
2£ + Aj channel, after the cuts given in [7] for 100 fb~ 1 of integrated luminosity. 



Mass (GeV) 


600 


800 


1000 


8M (D— > Z/W jet) 


120 GeV 


140 GeV 


160 GeV 


S 


133 


18 


2.5 


B 


9 


3 


2.8 


significance 


44.3 


6.8 


1.7 



The sum of number of signal and background events from both W and Z involving cases with calculated statistical 
significance, corresponding to one LHC year at high luminosity is presented in Table [5] Numbers of events in the 
signal and background events were calculated in the mass window of ±2 x 8M around the respective mass peak in the 
signal events. Since this decay channel has the biggest branching ratio among all D quark pair production processes, 
despite the obvious challenge of multi-jet environment, the statistical significance seems sufficient to observe D quarks 
up to a mass of around 800 GeV. 



3 Results and Conclusions 

The analysis presented here has shown that the lightest of the quarks predicted by the E§ GUT models can be dis- 
covered in different decay channels if its mass is less than 1 TeV By combining multiple channels from this study 
and from the previous one fOl . the discovery potential becomes higher and consequently the required luminosity 
becomes lower. Table [6] shows the expected number of signal and background events from each channel and the total 
significance after one year of running at nominal LHC luminosity of 100 fb _1 per year. To combine the significance 
of the different channels, we have used a simple method. For each channel, after obtaining the probability of compat- 
ibility with a background hypothesis using Poisson or Gauss statistics depending on the available number of events, 
we calculated the total probability for background compatibility of all four considered channels. The total probability 
was converted to significance assuming a normalized Gauss distribution. The top plot in Fig. [13] shows the combined 
significance for 30 and 100 fb~' of integrated luminosity as a function of the D quark mass. The same figure, on the 
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Table 6: Signal and background events for different channels and different masses of the D-quark for 100 ftr 1 inte- 
grated luminosity. Combined total significance is also given as the last row. 



mo(GeV) 


600 


800 


1000 


4£ + 2j Signal 
background 
-In p 


16 
3.0 
21.47 


3.7 
1.3 
4.78 


0.74 
0.4 
1.44 


2t + 2j + # T Signal 
background 
-In p 


53 
12 
120 


19 
13 
15.81 


4 
5 
3.32 


3£ + 2j + # T signal 
background 
-In p 


97 
24.9 
191.4 


18.3 
9.0 
20.66 


1.8 
1.05 
1.69 


2£ + 4j signal 
background 
-In p 


133 
9 
983 


18 
3 

25.3 


2.5 
2.8 
2.44 


— Eln p 

Combined Significance (a) 


1315.9 
51.3 


66.5 
11.3 


8.9 
3.8 



bottom part, contains the required luminosity in fb to observe 3(7 and 5(7 signals also as a function of the D quark 
mass. The required luminosity was calculated assuming Poisson statistics and the total number of events were rounded 
down to the nearest integer. The presence of at least one event was required for any given channel to be included in 
the total significance calculation. 




1 I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

500 600 700 800 900 1000 

m D (GeV) 




m D (GeV) 



Figure 13: Upper plot: Combined significance as a function of the D quark mass for an integral luminosity of 30 
(dashed line) and 100 (solid line) fb -1 . Lower plot: Required luminosities for the 3(7 observation limit (dashed line) 
and the 5(7 discovery limit as a function of the D quark mass. 

In conclusion, we show that the combination of the studied channels allows the discovery of the lightest isosinglet 
quark predicted by the Zsg GUT model within one year of LHC running time if its mass is up to about 950 GeV 
with a significance of 5(7 or more. As pair production was considered, these results are independent of the mixing 
angle between the D quark and its SM counterparts provided that the lightest iso-singlet quark mixes mainly to the 
first or second SM families. Although this study is based on a fast simulation of the detector response which was 
not fully validated and there are uncertainties associated with the QCD scale, statistical errors etc, we believe that 
the conservative selection cuts and the simplicity of the reconstruction algorithms give reliability to the conclusions. 
Additionally, the verification of these results using the full simulation of the detector is underway as preparation to 
first data from the LHC |T6). 
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